Dipole-dipole interaction between a quantum dot and graphene nanodisk 
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We study theoretically the dipole-dipole interaction and energy transfer in a hybrid system con- 
sisting of a quantum dot and graphene nanodisk embedded in a nonlinear photonic crystal. In our 
model a probe laser field is applied to measure the energy transfer between the quantum dot and 
graphene nanodisk while a control field manipulates the energy transfer process. These fields create 
excitons in the quantum dot and surface plasmon polaritons in the graphene nanodisk which interact 
via the dipole-dipole interaction. Here the nonlinear photonic crystal acts as a tunable photonic 
reservoir for the quantum dot, and is used to control the energy transfer. We have found that the 
spectrum of power absorption in the quantum dot has two peaks due to the creation of two dressed 
excitons in the presence of the dipole-dipole interaction. The energy transfer rate spectrum of the 
graphene nanodisk also has two peaks due to the absorption of these two dressed excitons. Addi- 
tionally, energy transfer between the quantum dot and the graphene nanodisk can be switched on 
and off by applying a pump laser to the photonic crystal or by adjusting the strength of the dipole- 
dipole interaction. We show that the intensity and frequencies of the peaks in the energy transfer 
rate spectra can be modified by changing the number of graphene monolayers in the nanodisk or 
the separation between the quantum dot and graphene. Our results agree with existing experiments 
on a qualitative basis. The principle of our system can be employed to fabricate nano-biosensors, 
optical nano-switches, and energy transfer devices. 

PACS numbers: 78.67.Hc, 73.20.Mf, 78.67. Wj 
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I. INTRODUCTION 



There has been growing interest in developing 
nanoscale optoelectronic devices by combining nano- 
materials with complementary optical properties into 
composite (hybrid) structures. The number of pos- 
sible composite systems that can be built from al- 
ready existing nanostructures is simply enormous. A 
significant amount of research on nanocomposites has 
been devoted to the study of exciton-plasmon inter- 
actions in metal-semiconductor nanostructures, which 
offer a wide range of opportunities to control light- 
matter interactions and electromagnetic energy flows on 
nanometer length scales. 1-6 Strong exciton-surface plas- 
mon coupling in semiconductor quantum dot (QD)-metal 
nanoparticle systems could lead to efficient transmis- 
sion of quantum information between qubits for applica- 
tions in quantum computing and communication. 2 These 
nanostructures also have applications in biophotonics 
and sensing, where nonradiative energy transfer between 
a QD and metal nanoparticle can be used to detect bio- 
logical molecules. 3 

In this paper we study theoretically the dipole-dipole 
interaction (DDI) and energy transfer between a quan- 
tum emitter and a graphene nanodisk. The quantum 
emitter can be a quantum dot, nanocrystal or a chem- 
ical or biological molecule. Here the quantum emittcr- 
graphene system is embedded in a photonic crystal, 
which acts as a tunable photonic reservoir for the emitter. 
Photonic crystals are engineered, periodically ordered 
microstructures that facilitate the trapping and control 
of light on the microscopic level. Applications for pho- 



tonic crystals include all-optical microchips for optical 
information processing, optical communication networks, 
sensors and solar energy harvesting. In our investiga- 
tion we consider a nonlinear photonic crystal, which has a 
refractive index distribution that can be tuned optically. 
The nonlinear photonic crystal surrounds the hybrid sys- 
tem and is used to manipulate the interaction between 
the quantum emitter and graphene nanodisk. 

Surface plasmon polaritons are created in the graphene 
nanodisk due to the collective oscillations of conduc- 
tion band electrons. They arise due to the dielec- 
tric contrast between graphene and the surrounding 
dielectric medium. Plasmonics is widely studied due 
to applications in ultrasensitive optical biosensing, 13 
photonic metamaterials, 14 light harvesting, 15 optical 
nanoantennas 16 and quantum information processing. 17 
Generally, noble metals are considered as the best 
available materials for the study of surface plasmon 
polaritons. 18 However, noble metals are hardly tunable 
and exhibit large Ohmic losses which limit their applica- 
bility to optical processing devices. Graphene plasmons 
provide an attractive alternative to noble-metal plas- 
mons, as they exhibit much tighter confinement and rel- 
atively long propagation distances. Furthermore, surface 
plasmons in graphene have the advantage of being highly 
tunable via electrostatic gating. Compared to noble met- 
als, graphene also has superior electronic and mechanical 
properties, which originate in part from its charge carri- 
ers of zero effective mass. 19 For example, charge carriers 
in graphene can travel for micrometers without scattering 
at room temperature. Graphene has also been recognized 
as a useful optical material for novel photonic and opto- 
electronic applications. 20-25 For these reasons, the study 
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of plasmonics in graphene has received significant atten- 
tion both experimentally and theoretically. 21 ' 22 ' 26-29 

Recently, experimental research on graphene has been 
extended to the fabrication and study of QD-graphene 
nanostructures. 30-34 For example, a CdS QD-graphene 
hybrid system has been synthesized by Cao et al., 30 
in which a picosecond ultrafast electron transfer pro- 
cess from the excited QD to the graphene matrix was 
observed using time-resolved fluorescence spectroscopy. 
Chen et al. 31 have fabricated CdSe/ZnS QDs in contact 
with single- and few-layer graphene sheets. By measur- 
ing the strong quenching of the QD fluorescence, they 
determined the rate of energy transfer from the QD to 
graphene. A similar study by Dong et al. 32 was per- 
formed on a CdTe QD and graphene oxide system, but 
in their case the QDs were modified with molecular bea- 
cons in order to demonstrate that the hybrid system 
can be used for sensing biological molecules. Wang et 
al. 33 have synthesized graphene-CdS and graphene-ZnS 
QD hybrid systems directly from graphene oxide, with 
CdS and ZnS QDs very well dispersed on the graphene 
nanosheets. They also measured the QD photolumines- 
cence and observed the energy transfer between the QDs 
and graphene. Metal nanoparticle-graphene hybrid sys- 
tems have also been fabricated by several groups. 21 ' 35 " 37 

Here we study a QD-graphene hybrid system, where 
energy transfer occurs due to the interaction between 
optical excitations in the QD and graphene nanodisk. 
The optical excitations in the QD are excitons, which 
are electron- hole pairs, while those in the graphene nan- 
odisk are surface plasmon polaritons, which are created 
due to the collective oscillations of conduction band elec- 
trons. The QD is taken as a three-level lambda-type 
system in which two distinct excitonic transitions occur. 
Other three-level systems in the ladder- 3 and V-type 38 
configurations interacting with a metallic nanoparticle in 
the presence of two external fields have been studied. In 
our model we include a probe laser field which is coupled 
with one excitonic transition and measures the energy 
transfer spectra of the QD and graphene. We also con- 
sider that a control laser field is applied to monitor and 
control the energy transfer. Aside from creating excitons 
in the QD, these fields also generate surface plasmon po- 
laritons in graphene. The dipoles created by excitons in 
the QD and plasmons in the graphene nanodisk then in- 
teract via the DDI. This interaction is strong when the 
QD and graphene are in close proximity and their optical 
excitation frequencies are resonant. 

We have found that the power absorption spectrum of 
the QD has two peaks when the QD and graphene nan- 
odisk are in close proximity, indicating the creation of 
two dressed excitons due to the DDI. These dressed exci- 
tons are transported to graphene, and produce two peaks 
in the spectrum of the energy transfer rate to graphene. 
We show that the energy transfer between the QD and 
graphene can be switched on and off by changing the 
strength of the DDI coupling or by applying an intense 
laser field to the nonlinear photonic crystal. The inten- 



sities of peaks in the energy transfer rate spectra can be 
controlled by changing the number of graphene mono- 
layers or by changing the distance between the QD and 
graphene. We have also predicted that the intensity of 
these peaks can be modified in the presence of biolog- 
ical materials. Our findings agree with the experimen- 
tal results of Refs. 30-34 on a qualitative basis. The 
present system can be used to fabricate nano-biosensors, 
all-optical nano-switches and energy transfer devices. 



II. THEORETICAL FORMALISM 

We investigate theoretically the dipole-dipole interac- 
tion and energy transfer between a quantum dot and 
graphene nanodisk when the system is embedded in a 
nonlinear photonic crystal. A schematic diagram for the 
present system is shown in Fig. 1. A graphene nanodisk 
(or nanoflake) lies in the x-y plane, on top of which a QD 
is deposited. The center-to-center distance between the 
QD and graphene nanodisk is taken as R. The combined 
QD-graphene nanodisk system can also be referred to as a 
QD-graphene nanocomposite or hybrid system. The QD 
considered here has three discrete states, where 2} and 
1 3) are the lower-energy states which are near-degenerate 
and both coupled to the common optically excited state 

I). The QD acts as a three-level quantum emitter, where 
excitons are created by the transitions |2) o |1) and 

3) <H> 1 1) with resonance frequencies (dipole moments) 
W12 (M12) an d w i3 (^13)1 respectively. This so-called 
Lambda-type energy level configuration has been widely 
studied in atoms, where quantum optical effects such as 
electromagnetically induced transparency and coherent 
population trapping have been demonstrated. 39 ' 40 More 
recently, this energy level configuration has been achieved 
in semiconductor QDs. 41-43 

In our model, we consider that a probe field E2 = 
-E°cos(w2£) is applied between states |1) and |2). To 
control the coupling between the QD and graphene a 
control field E 3 = E® cos (w 3 i) is applied between states 

I) and |3). These laser fields excite both the QD and 
graphene nanodisk. In the QD, the probe and control 
fields create excitons which produce dipole electric fields 
that interact with the graphene nanodisk. Similarly in 
the graphene nanodisk surface plasmons are generated 
by the probe and control fields, which produce dipole 
electric fields that interact with the QD. 

The dipole electric fields produced by the QD and 
graphene are written as 
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respectively. Here thd = (2efc + e<i) /3efc, where e<j is the 
dielectric constant of the QD and e b is the dielectric con- 
stant of the medium surrounding the QD-graphene sys- 
tem. The parameter gi (I — x, y, z) is the polarization 
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where ai is the polarizability of the graphenc nanodisk 
and e g (u>) denotes the dielectric function of graphene. 
Here Tz is called the depolarization factor, and is obtained 
as 



?! = 



L x LyL z 



dq 



(Lf + q)f(q) 



(6) 



where 

/ (?) = y/(Ll + q) (Ll + q) (Ll + q) (7) 

The depolarization parameters satisfy the relation q x + 
<; y + <; z = 1, and determine the optical response of the 
graphene nanodisk based on its shape. 44 For an oblate 
spheroid where L x — L y and L x > L z , the depolarization 
factors reduce to 
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FIG. 1: A schematic diagram of the QD-graphene nanocom- 
posite embedded in a photonic crystal. The QD has three dis- 
crete states, where |2) and |3) denote the lower-energy states, 
which are near-degenerate and are both coupled to the com- 
mon optically excited state |1). 



parameter, with g x = g y = — 1 and g z = 2 for electric 
fields polarized in the x-y plane or in the z-direction, 
respectively. In Eq. (1) Pqd and Pq are the polariza- 
tion of the quantum dot and graphene nanodisk, respec- 
tively. The polarization of the QD is obtained as 



Pqd = Vii (Pu+Pn) 



(2) 



i=2,3 



where ^ and are the dipole moment and density 
matrix element, respectively, for the transition \i) o \j). 

The total electric field felt by the graphene nanodisk 
is written as 



Eq — E 2 ~t~ E3 -\- E 



QD 
DDI 



(3) 



where the first, second and third terms are the contribu- 
tions from the probe, control and the QD dipole field, re- 
spectively. Using the quasistatic dipole approximation, 44 
the polarization in the graphene nanodisk is 



QD 
DDI 



P G = e b ai[E 2 +E z + E 

AttL x L v Lz [e m (u>) - e b 



at = 



3e b + 3<^ [e m (w) - e b ] 



(4) 
(5) 



where the eccentricity of the nanodisk is defined as 



1 — (L x /L z ) ■ Note that this equation is valid for both 
cases where L x > L z and L x < L z - A4l 

The above expression for the polarizability given in Eq. 
(5) has been widely used in the literature to study the op- 
tical properties of metallic nanodisks, and has been found 
to give good agreement with experimental results. 45,46 
In the quasistatic approximation the dimensions of the 
graphene nanodisk are much smaller than the wavelength 
of incident light and we therefore assume a spatially uni- 
form but time-varying electric field across the graphene 
nanodisk. Here the wavelengths of light considered are 
on the order of several hundred nanometers, and thus the 
size of the graphene nanodisk must be less than 100 nm. 
It is important to note that our model is only valid for 
nano-size graphene samples and not for bulk materials. 
For a very flat and thin disk we take L x ^> L z , which 



gives 



47,48 
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The above method has been used for a graphene flake in 
Rcf. [49]. 

The total electric field felt by the QD is written as 



E, 



QD 



E 2 E 3 9iai [E 2 + E 3 + E% u Dl 
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where the first, second and third terms are the contri- 
butions from the probe, control and the graphene dipole 
field, respectively. Putting the expression of E c j^ )I from 
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Eq. (1) into the above expression we get 



Eqd — 



^ ( E F + E DDl) 

i=2,3 
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The dipoles of the QD interact with the dipole elec- 
tric field produced by the graphene nanodisk, and vice 
versa. This interaction is called the dipole-dipole inter- 
action (DDI). The terms in the interaction Hamiltonian 
of the QD due to the external probe and control fields 
and the DDI are expressed in the rotating wave approx- 
imation as 



Hqd-f = - J2 VuEUti + fc-c (13) 

i=2,3 

Hqd-ddi = - ^ Li u E DDI a++ h.c. 



i=2,3 



Using the expressions for E F and E DDI in Eqs. (11) 
in the above expression and putting the Hamiltonian in 
interaction representation we get 

H Q d-f = - ^«4e" i(ui " Wli)< + h - c - ( i4 ) 

i=2,3 

Hqd-ddi = -^ft(n,+ Aip^a+e-^-^ + h.c. 

i=2,3 

where cr^- = \i) (j\ (<Jij = \j) (i\) is the exciton cre- 
ation (annihilation) operator. The interaction Hamilto- 
nian term Hqd-f given in Eq. (14) represents the direct 
contribution from the external probe (i ~ 2) and control 
(i = 3) fields incident on the QD, and is given in terms 
of the Rabi frequencies tti which measure the intensities 
of these fields. The second contribution Hqd-ddi rep- 
resents the fields incident on the QD due to the DDI be- 
tween the QD and graphene nanodisk, and contains two 
terms Ui and Aip u . The first DDI term is due to the 
interaction of the QD with the dipole electric fields from 
the graphene nanodisk induced by the probe and con- 
trol fields. Therefore we refer to this as the direct DDI 
term. The second DDI contribution is due to the inter- 
action of the QD with a dipole field from graphene that 
arises when the external fields polarize the QD, which 



in turn polarize graphene. In other words these contri- 
butions are the self-interaction of the QD, as they de- 
pend on the polarization of the QD. For this reason this 
term is called the self-induced DDI parameter. The sur- 
face plasmon polariton resonance frequency uj l sp in the 
graphene nanodisk is obtained by setting the real part of 
the denominator in ai (oj) equal to zero and solving for 
u. When the optical excitation frequencies of the QD lie 
near the surface plasmon polariton resonance frequencies 
of the graphene nanodisk (i.e. lou « w sp ), the DDI be- 
comes very strong due to the enhanced local fields in the 
vicinity of the graphene nanodisk. This interaction leads 
to excitation and energy transfer between the QD and 
graphene. 

The combined QD-graphene system is embedded in a 
photonic crystal consisting of dielectric spheres arranged 
periodically in three dimensions, which acts as a reser- 
voir for the QD. Therefore, we consider that the excited 
state |1) spontaneously decays to the lower-energy states 
|2) and |3) due to excitons coupling with Bloch photons 
in the photonic crystal (see Fig. 1). In this case the 
interaction Hamiltonian is given as: 

Hqd-pc = - Y.^^e^-^ + h.c, (15) 



9u 



i=2,3 k 



n^k 

2e b V PC 
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where is the polarization unit vector of the Bloch pho- 
tons in the PC and Vpc is the volume of the photonic 
crystal. The operator (dk) is the photon creation (an- 
nihilation) operator, while ujk and k are the Bloch photon 
frequency and wave vector, respectively. 

Finally, the total interaction Hamiltonian of the system 
is written as 



H 



Hqd-f + Hqd-ddi + Hqd-pc (16) 



i=2,3 



i=2,3 fc 



where 



9i = arctan 



Re(n,)] 2 + [im(ni)] 2 



ImOT) 



fii +ReOT) 



Note that the first term in the above expression for H is 
the combination of Hqd-f with the direct DDI term. 

We use the density matrix method to evaluate the en- 
ergy transfer between the QD and the graphene. Using 
Eq. (16) for the interaction Hamiltonian and the master 
equation for the density matrix 39 we obtained the fol- 
lowing equations of motion for the QD density matrix 
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elements as 
dp 22 



dt 



= 2r 2 ip n - iR 2 e l 2 p 21 - iA 2 p 12 p 21 



+ iR 2 e- t82 p 12 + iA. 
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1 - d i3Pi3 + iR2e i02 P 23 + iA 2Pi 2 p 23 
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iR 3 e l03 p 21 - iA 3 p 13 p 21 



du — r 2 i + r 3 i 

Si = OJi — LO\i. 



id 
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Here <5j are the detuning of the probe (i = 2) and control 
(i = 3) fields. Note that the diagonal elements of the 
density matrix satisfy the relation p ±1 + p 22 + p 33 = 1. 
The quantities and are the non-radiative decay 
rate and energy shift, respectively due to self-induced 
DDI parameters Aj. They are found as 

T id = Im(Aj) {p u - p n ) 
A id = Re(A;) - p n ) 

The parameters Tn represent the spontaneous decay 
rates of excited state |1) to state |i) due to the Bloch 
photons in the photonic crystal, and are given as 



-D{u u ) 



where 
and 
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In the above expression, is the exciton decay rate 
due to the background radiation field in free space. Here 
L = 2a + 2b is the photonic crystal lattice constant, 2b 
is the spacing between dielectric spheres and a is the ra- 
dius of the spheres. Parameters n a and n b denote the 



refractive index of the dielectric spheres and background 
material in the photonic crystal, respectively. The ex- 
pression for the photonic density of states D(u>) has been 
derived in Ref. [50]. Here we have used the Marko- 
vian approximation in order to derive the decay rates for 
the QD in the presence of the photonic crystal. This 
approximation ignores memory effects in the electromag- 
netic reservoir due to the presence of the QD, and is 
valid when the photonic density of states can be consid- 
ered smooth and slowly- varying compared to the energy 
difference between the edge of the photonic band gap 
and the resonance frequency of the QD. 51 Note that in 
our calculations, we remain within the regime where the 
Markovian approximation is valid. Therefore the effect of 
the photonic crystal serves only to alter the decay rates of 
the excitonic transitions compared to those in free space. 
We numerically solve Eqs. (17-21) by first substituting 
P12 = 7>i2^\ Px3=~Px3^ and p 23 = ~p 23 e-^- e -\ 

Following the method of Ref. [44] and using Eq. (16), 
the energy absorption rate of the QD (Wqd) and the 
energy transfer rate from the QD to graphene (Wq) are 
found as 



Wqd = ^ fruup^rn 



i=2,3 



w G =j2 



9l, z p\iUi Im(a x ) \p u 



tts 87r2e 6 e 2 Je 6g |^6 



(23) 
(24) 



where e bg = (2e b + e g ) /3e b . The expression for Wqd is 
obtained by assuming that the power radiated from the 
QD is equal to its energy absorption rate. Similar expres- 
sions have been widely used in the literature on hybrid 
systems. 2,3 Note that the energy transfer to graphene de- 
pends on the coherences p u of the QD density matrix, 
which change depending on the center-to-center distance 
between the QD and graphene, R. Therefore Eq. (24) 
does not simply depend on R 6 , but rather is a much more 
complicated function of R. 



III. RESULTS AND DISCUSSION 

We consider a graphene nanodisk with thickness L z = 
0.35 nm (i.e. a single graphene layer) and size ratio 
L x /L z = 20. The plasmon frequency and background 
dielectric constant of graphene arc taken from experi- 
ments as 6.02 eV and 1.964, respectively. 21 The decay 
rate in graphene is taken as j G = 5 THz, which is con- 
sistent with the relaxation rates reported in refs. [22] 
and [26]. With these parameters, the surface plasmon 
polariton resonance frequencies in the graphene nanodisk 
are calculated as Tiu: x sp — 0.8026 cV and huj z sp = 4.1250 
cV. The QD dielectric constant and dipole moments are 
taken as e d = 12 and /i 12 = p 13 = 0.1 e x nm, re- 
spectively, while the free space decay rates for the QD 
are taken as r 21 = T 31 = 0.2 peV. These parameters are 
comparable to those commonly found in the literature for 
QDs. 2 ~ 4 Here the transition energies in the QD are taken 
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to lie near the plasmon resonance htu x p as Tiuj\2 — 0.8046 
cV and hu>is = 0.8036 eV. The combined QD-graphene 
nanodisk hybrid is contained within a photonic crystal 
made of polystyrene spheres arranged periodically in air. 
Similar photonic crystals have been fabricated by Liu et 
al. 9 , in which ultrafast all-optical switching was experi- 
mentally demonstrated. Photonic crystal parameters are 
taken as a — 170 nm, L = 480 nm, n a — 1.59 and n& = 1. 
With these parameters, we find that a photonic band gap 
appears between frequencies 0.8225 eV and 0.9843 eV. 
Note that the lower edge of the band gap lies near oj x p and 
the QD transition frequencies W12 and W13. The vacuum 
decay rates for the QD are taken as r° = r§ = 0.2 /xeV, 
and in the presence of the photonic crystal it is found 
that T 2 i = 1.1370 /icV and T 3 i = 1.1127 fj,eV. Here the 
background dielectric constant was taken as e# = 2.081. 
Throughout the following calculations we consider that 
the intensity of the probe and control fields are 1.0 and 
3.0 W/cm 2 , respectively. 

We first consider the case where the excitonic transi- 
tion |3) <-> |1) is coupled with the surface plasmon res- 
onance of the graphene nanodisk. In this configuration, 
the transition frequency a; 13 is near uj x p while both the 
control field E3 and the transition dipolc moment /i 13 
are polarized in the x-y plane. Conversely, the transition 
|2) "R- |1) is not coupled with the graphene nanodisk. 
This situation occurs when the probe field E 2 and tran- 
sition dipole moment fi 12 are polarized in the z-direction 
and CJ12 is far away from Lo z sp (see Fig. 2 inset). The 
energy absorption rate in the QD is evaluated from Eq. 
(23) and the results are presented in Fig. 2a when the 
QD-graphene separation R is varied and the control field 
is resonant with the |3) <-» |1) transition such that <5 3 = 0. 
It is found that the power absorption spectrum has a sin- 
gle peak with an extremely narrow transparent window 
at S 2 = when the QD and graphene are further away 
with each other (i.e. R — 20 nm). This narrow minima is 
due to electromagnetically induced transparency in the 
system. When the QD is close to graphene (i.e. R = 8 
nm) the power absorption peak splits into two peaks and 
a clear minima appears at S2 = 0. The observed split- 
ting is due to the DDI and surface plasmon polariton 
coupling. 

The splitting of the power absorption spectrum can 
be explained using the theory of dressed states. When 
the QD is close to the graphene nanodisk there is strong 
coupling due to the DDI for the transition |3) <H> |1). 
This causes the excited state |1) to split into two dressed 
states, namely |1 + ) and |1-}. Therefore, there arc now 
two transitions |2) <-» |1 + ) and |2) o |1_) which give two 
peaks and and a minima in the spectrum. In other words, 
a single exciton splits into two dressed excitons, and their 
energy difference is found to be proportional to the DDI. 
As the distance between the QD and graphene increases, 
the splitting decreases since the direct DDI term IT is 
inversely proportional to i? 3 . In Fig. 2b we have plotted 
the energy absorption rate when the control field is de- 
tuned such that S3 = 10 /xeV. In this case, we find one 



FIG. 2: Energy absorption rate of the QD as a function of 
probe field detuning 82 when the QD-graphene nanodisk sep- 
aration R is varied, (a) £3 = 0; (b) S3 — 10 fieV. Inset: 
Polarization of the probe and control fields. 



peak and negligible electromagnctically-induced trans- 
parency when R is large. When R decreases, the single 
peak splits into two peaks due to the DDI. These results 
show that the DDI can be used to split one exciton into 
two excitons, and also to control the electromagnetically- 
induccd transparency phenomenon. 

In Fig. 3 we have investigated the effect of the photonic 
crystal on the energy absorption rate in the lambda-type 
QD. Initially, the lower band edge of the photonic crys- 
tal lies far away from the resonance energies of the QD 
(see solid curve), and there is weak coupling between the 
QD and photonic crystal. When we move the lower band 
edge closer to the resonance frequency W13 of the QD, the 
two peaks in the power absorption spectrum merge into 
a broad peak with a narrow electromagnetically-induccd 
transparency window at S2 = (see dashed curve). Note 
also that the height of the peaks decreases. The merg- 
ing of the split peaks in the QD power absorption spec- 
trum occurs because the spontaneous decay rates become 
larger than the DDI splitting for the two peaks. The 
value of the decay rate is large because the photonic den- 
sity of states is large when the resonance energy of the 
QD lies near the band edges. For example, we found 
T 2 i = 6.40 fieV and T 3 i = 3.81 fieV whereas the energy 
splitting is about 2.80 /xeV. Here, the location of the pho- 
tonic crystal band edges can be changed by applying an 
intense pulsed laser field. The intense laser field causes 
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FIG. 4: Energy transfer rate from the QD to graphene as a 
function of probe field detuning 82 when the QD-graphene 
nanodisk separation R is varied and £3 = 0. Inset: Polariza- 
tion of the probe and control fields. 



FIG. 3: Energy absorption rate in the QD as a function of 
probe field detuning 82 when the lower band edge of the pho- 
tonic crystal is taken as e v = 7iwi2 + 17.88 meV (a) and 
e v = %uj\2 + 0.56 meV (b). Here R — 13 nm and S3 = 0. 
Inset: Polarization of the probe and control fields. 



the refractive index of polystyrene, an optical nonlinear 
material, to change due to the Kerr effect. This change 
is quantified by the expression n' a = n a + n^I pump where 
ri3 is the Kerr nonlinearity constant and has the value 
n 3 = 1.15 x 10~ 12 cm 2 /W for polystyrene. 10 For the 
pump field intensity I pump = 31.0 GW/cm 2 we found 
that the photonic crystal band edge shifts such that 
Ae^ = — 17.32 meV. This means that the hybrid system 
can be used to study the nonlinear properties of photonic 
crystals. Using an external pump field to induce a large 
Kerr nonlinearity in the polystyrene photonic crystal is 
also an effective way to switch the energy transfer be- 
tween two states; from high to low energy transfer peaks. 
Alternatively, the refractive index of the background ma- 
terial in the photonic crystal can also be modified by im- 
mersing the photonic crystal in another material. There- 
fore the present QD-graphene system can also be used as 
a nano-sensor. 

We now investigate the energy transfer rate from 
the QD to graphene, and the results are plotted in 
Fig. 4 as a function of probe detuning when the QD- 
graphene separation R is varied. We find that the 
power transfer spectrum has a single peak with a narrow 
electromagnetically-induced transparency window when 
R is large (i.e. R = 20 nm). When the QD is brought 
closer to graphene (i.e. R = 8 nm), the power trans- 
fer spectrum has one large minima and two peaks with 
separation proportional to the DDI. This indicates that 
energy is transferred from the QD to the graphene when 



the two dressed excitons created in the QD are absorbed 
by graphene. This is an interesting finding, and can be 
used to transfer energy absorbed by the QD from a light 
source (i.e. the sun) to graphene where it can be stored. 
Therefore, one can fabricate energy transfer and storage 
devices (i.e. solar cells) from the present hybrid system. 

It is also found that the height of the energy transfer 
peaks increases as the QD-graphene separation decreases 
(see Fig. 4). This effect has been observed experimen- 
tally by Chen et al. 31 and Dong et al. 32 They found that 
as the distance between CdTe-QDs and a graphene oxide 
sheet decreases there is a strong quenching of the QD flu- 
orescence. They concluded that the fluorescence quench- 
ing could be due the energy transfer from the QD to the 
graphene sheet. For example, Chen et al. 31 deposited 
graphene on quartz substrates and then CdSe/ZnS-QDs 
were deposited on graphene. The fluorescence measure- 
ments were performed on the individual QDs located 
both on the bare quartz substrate and on a graphene 
layer. They observed strong fluorescence quenching for 
QDs deposited on the graphene sheet, which was at- 
tributed to the energy transfer between QD and graphene 
and not due to photoinduced electron transfer from the 
QD to graphene. Similarly Wang et al. 33 performed pho- 
tolumincscence measurements on CdS-QDs and ZnS 
QDs on graphene and observed a strong quenching of 
photoluminescence for these QDs due to the presence of 
the graphene sheet. They also performed transient pho- 
tovoltaic experiments on their hybrid systems and found 
a very unexpected strong positive photovoltaic response 
due to the DDI. Conversely, it was found that separate 
samples of graphene and CdS-QDs of a similar size do not 
show any photovoltaic response. They concluded their ex- 
perimental findings can be explained due to the energy 
transfer between the QD and the graphene sheet. Similar 
energy transfer between a QD and carbon nanotube has 
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also been found experimentally by Shafran et al.. 34 

When a QD is in contact with biomolecules, molecular 
beacons, DNA or aptamers, its dielectric constant can be 
modified. Therefore we have investigated the role of the 
dielectric constant of the QD on the energy transfer to 
graphene. The results are plotted in Fig. 5a for three val- 
ues of td ■ It is found that by changing the dielectric con- 
stant of the QD, the height of the energy transfer peaks 
can be modified. For example by increasing or decreasing 
the dielectric constant the height of the energy transfer 
spectra decreases or increases, respectively. This is be- 
cause the the energy transfer is inversely proportional 
to the square of the dielectric constant, as shown in Eq. 
(24). This effect has also been verified experimentally by 
Dong et al., 32 where upon integrating a molecular bea- 
con to a CdTe-QD it was found that the fluorescence 
quenching due to graphene is modified. We also note 
that at certain values of probe detuning, say for example 
82 ~ ±1.5 fieV, the sensitivity of the energy transfer rate 
to the change in dielectric constant is quite high. This is 
an interesting finding, particularly if one considers that 
the present hybrid system can be used to fabricate nano- 
biosensors. 

We have also studied the effect of the number of 
graphene layers on the energy transfer spectrum. In Fig. 
5b we have plotted the energy transfer rate to graphene 
when a single graphene layer or two layers are considered. 
Here the ratio L x /L z = 20 is preserved in order to keep 
the surface plasmon polariton resonance frequency con- 
stant. Note that for two layers of graphene the height 
of the energy transfer peak increases. As we add ad- 
ditional layers of graphene, we increase its volume. In 
turn, the DDI between the QD and graphene is enhanced. 
Therefore, both the height of the peaks in the energy 
transfer spectrum and their splitting are increased. This 
effect has also been verified experimentally by Chen et 
al. , where it was found that increasing the number 
of graphene layers in a CdSc/ZnS nanocrystal-graphene 
composite system enhanced the QD fluorescence quench- 
ing effect. 

We now consider an alternative configuration for the 
QD-graphene nanocomposite system where both transi- 
tions |2) O |1) and |3) O |1) couple with the surface 
plasmons in the graphene nanodisk. In this configura- 
tion we consider that both W12 and w 13 are close to uo x sp 
and the transition dipole moments (fields) /i 12 (probe 
field) and /i 13 (control field) are aligned along the x- and 
y-directions, respectively. In Fig. 6a we have plotted the 
energy transfer rate from the QD to graphene as a func- 
tion of the probe field detuning while varying R. Here 
the physical parameters are the same as considered in 
our previous calculations. Note that we see two peaks 
due to the DDI effect as in the first configuration. Pre- 
viously the two peaks were symmetric, but in this case 
they are asymmetric. This is due to the self-induced DDI 
parameter A 2 , which causes both peaks to shift towards 
positive detuning due to the change in the effective probe 
detuning from 62 to 62+ A2<j, as shown in Eq. (1191) . Here 



FIG. 5: (a) Energy transfer rate from the QD to graphene 
when the dielectric function of the QD is taken as = 10 
(dotted curve), 12 (solid curve) and 14 (dashed curve). Here 
R — 13 nm and £3 = 0. (b) Energy transfer rate from the QD 
to graphene when the thickness of graphene is varied between 
one layer (solid curve) or two (dashed curve). Here R = 13 
nm and S3 = 0. Inset: Polarization of the probe and control 
fields. 



also the width of both peaks increases due to the non- 
radiative decay T2d- In the previous configuration the 
self-induced DDI parameter A 2 was zero because there 
was no coupling between the QD transition 2} <->■ |1) 
and graphene. These effects are enhanced by decreasing 
R. 

In Fig. 6b the effect of the photonic crystal is investi- 
gated in the same way as in Fig. 3, and similar results 
are found as for the previous configuration. By applying 
an external pump laser field to the polystyrene photonic 
crystal, the power transfer to graphene can be switched 
from high to low values. We note that due to the asym- 
metry of the power transfer spectrum in this configura- 
tion, the sensitivity of this switching effect can change 
drastically depending on the value of probe field detun- 
ing. For example, negative detunings close to S2 = 
show a sharp peak in the energy transfer spectrum when 
the pump field is absent, and this peak is suppressed 
when the pump field is applied. 
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FIG. 7: Energy absorption rate of to the ladder-type QD as 
a function of probe field detuning 82 when the QD-graphene 
nanodisk separation R is varied. Here hu)23 = 0.80 36 eV 
and the intensities of the probe and control fields are 1.0 and 
3.0 W/cm 2 , respectively. Other parameters are the same as 
considered previously. Inset: Schematic of the QD-graphene 
hybrid system with ladder-type energy level structure. Here 
DDI coupling occurs only for the |2) o |3) transition. 



FIG. 6: Energy transfer rate from the QD to graphene as a 
function of probe detuning S2 for the second configuration of 
dipole moments and fields (see inset), (a) The QD-graphene 
separation is varied from R = 13 nm (solid curve) to R = 11 
nm (dashed curve), (b) R = 11 nm and the lower band edge 
of the photonic crystal is taken as e v = hu)\2 + 17.88 meV 
(solid curve) and e v = ?koi2 + 0.56 meV (dashed curve). Here 
83 = 0. Inset: Polarization of the probe and control fields. 



feet of the photonic crystal on the energy absorption rate 
in the ladder-type QD, and the results are shown in Fig. 
8. Again we consider that a pump field of intensity 31.0 
GW/cm 2 is applied, which causes the photonic crystal 
band edge to shift and increases the decay rate of the QD. 
Note that the power absorption spectrum merges into a 
broad peak in the same way as it did for the lambda-type 
QD (see Fig. 3b). Here, however, the narrow minima 
present in the lambda-type QD is absent and we clearly 
see two peaks merging into one. 



IV. CONCLUSIONS 



Finally, we investigate the energy absorption rate in 
a ladder-type QD coupled with the graphene nanodisk. 
The formulation for this system is given in the appendix. 
We consider the case where the control field is coupled 
with the QD transition 2} «-> |3) and the graphene nan- 
odisk, while the probe field is only coupled to the QD 
transition |1) «-> 2). Here the resonance frequency U23 
lies near wj p , while UJ12 is uncoupled from both Wg p and 
u z sp . This situation is analogous to that explored in Fig. 
2 for the lambda-type QD. In Fig. 7 the energy absorp- 
tion rate in the ladder-type QD is plotted when the QD- 
graphene nanodisk separation R is varied. We find that 
the power absorption spectrum gives two peaks and a 
minimum when R is small (i.e. R = 8 nm), as was found 
in Fig. 2. Note that for the ladder- type QD, the narrow 
minima due to electromagnetically induced transparency 
does not appear. This is because the same electromag- 
netically induced transparency effect does not appear in 
ladder-type systems. 39 We have also investigated the ef- 



We have investigated the dipolc-dipole interaction 
(DDI) and energy transfer in a quantum dot (QD)- 
graphene nanodisk system embedded within a photonic 
crystal. Our simulations predict that in this system, mul- 
tiple excitonic states (dressed states) can be created in 
the quantum dot and then transferred to graphene with 
different frequencies. This phenomenon occurs purely 
due to the DDI between the QD and graphene, and re- 
sults in energy transfer. We have demonstrated that the 
energy absorption of the QD and/or the energy transfer 
from the QD to graphene can be switched on and off by 
changing the strength of the DDI or by applying an in- 
tense external field to the photonic crystal. We have also 
verified our findings qualitatively with recent experimen- 
tal data on energy transfer in QD-graphene nanocompos- 
ite systems. Our numerical results provide motivation 
for future experimental and theoretical investigations on 
nanocomposites made from graphene, carbon nanotubes, 
quantum dots and photonic crystals. The present theory 
can be applied to hybrid systems consisting of graphene 
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FIG. 8: Energy absorption rate of the ladder-type QD as 
a function of probe field detuning 82 when the lower band 
edge of the photonic crystal is taken as e v = 7iW23 + 17.88 
meV (a) and e v = huj23 + 0.10 meV (b). Here R = 10 nm, 
S3 = and other parameters are the same as in Fig. 7. Inset: 
Schematic of the QD-graphene hybrid system with ladder- 
type energy level structure. Here DDI coupling occurs only 
for the 1 2) O |3) transition. 



with quantum emitters such as quantum dots, nanocrys- 
tals, atoms and chemical or biological molecules; the only 
requirement is the quantum emitter should have at least 
three states. The proposed nanocomposite system can be 
used to fabricate nano-sensors, all optical nano-switches, 
energy transfer devices and energy storage devices. 



V. APPENDIX 

We consider a three-level quantum dot in the ladder 
configuration, where |1), 2) and |3) denote the ground, 
first excited and second excited states, respectively. The 
probe field with amplitude E 2 and frequency uj 2 is cou- 



pled between states |1) and |2), while the control field 
with amplitude E 3 and frequency uj 3 is coupled between 
states |2) and |3). The decay of level |2) to level |1) (level 
|3) to level |2)) is given as T 2 i (r 32 ). Using the same 
methods as for the lambda-type system, the density ma- 
trix equations of motion for the ladder-type energy level 
configuration are then obtained as 



dp 



22 



dt 



dp 



33 



dt 

dp 2 i 
dt 

d P 32 

dt 

dp 3 i 
dt 

where 



and 



-T 21 p 22 + r 32 p33 + iR 2 e 162 p 21 + tA* 2 p 12 ^l) 

-iR 2 e l02 p 12 - iA 2 p 21 p 12 - iR 3 e~ 103 p 32 
-« A aP23P32 + iR3e lf>3 p 23 + iA 3 p 32 p 23 

-r 3 2P 3 3 + ^3e^ 93 /?32 + «A*p 23 p 32 (A2) 
-iR 3 e l03 p 23 - ik 3 p 32 p 23 

d 2 ip 21 + iR 2 e l ° 2 (p 22 - p n ) - iR 3 e~ m p 31 (A3) 

-* A 3P23P31 

dz2P 32 + iR3e m (p 33 - p 22 ) + iR 2 e- M2 p 31 (A4) 

+« A 2Pl2P31 

(iS 2 + iS 3 - r 32 ) p 31 - iR 3 e' 93 p 21 - iA 3 p 3 £££) 
+iR 2 e l62 p 32 +iA 2 p 21 p 32 



d 2 \ 
d 32 



15 2 + iA 2d - 

15 3 + iA 3d 



r 2 i 



r 2 d 
r 3 d 



A 2d 

r 2 d 

A 3d 

F 3d 



Re(A 2 ) (p 22 -p n ) 
Im (A 2 ) (p 22 - p n ) 

R- C ( A 3)(P 3 3-P2 2 ) 

Im (A 3 ) (p 33 - p 22 ) 



In the above expressions, all quantities are the same as 
given previously for the lambda-type system but with the 
substitutions cji 3 — > lu 23 , p 13 — > p 23 and T 31 — > T 32 . The 
QD energy absorption rate and the power transfer in this 
system are calculated using Eqs. (23) and (24) with the 
subsitutions p n ->■ p u , p 13 p 23 . 
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